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Abstract 

The ultracold molecular conversion rate occurring in an adiabatic ramp through a Fano-Feshbach 
resonance is studied and compared in two statistical models. One model, the so-called stochastic 
phase space sampling (SPSS)[E.Hodby et at, PRL.94 120402(2005)] evaluates the overlap of two 
atomic distributions in phase space by sampling atomic pairs according to a phase-space crite- 
rion. The other model, the chemical equilibrium theory (ChET)[S.Watabe and T.Nikuni, PRA.77 
013616(2008)] considers atomic and molecular distributions in the limit of the chemical and ther- 
mal equilibrium. The present study applies SPSS and ChET to a prototypical system of K+K— > 
K2 in all the symmetry combinations, namely Fermi-Fermi, Bose-Bose, and Bose-Fermi cases. To 
examine implications of the phase-space criterion for SPSS, the behavior of molecular conversion is 
analyzed using four distinct geometrical constraints. Our comparison of the results of SPSS with 
those of ChET shows that while they appear similar in most situations, the two models give rise 
to rather dissimilar behaviors when the presence of a Bose-Einstein condensate (BEC) strongly 
affects the molecule formation. 

PACS numbers: 31.15.bt,67.85.Pq,82.60.Hc 
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I. INTRODUCTION 



In recent years, the production and application of ultracold molecules have attracted 
much attention. Ultracold polar molecules in particular are being investigated eagerly be- 
cause of their great promise for numerous applications such as the quantum simulator l|, 



ultra high precision measurements |2|], ultracold chemical react ions |3|], etc. However, adap- 
tation of standard laser cooling techniques to molecules is not readily achieved. Instead 
experimentalists have employed a scheme to convert ultracold atoms into molecules without 
heating, i.e. the combination of Fano-Feshbach resonance (FFR) and Stimulated Raman Adi- 
abatic Passage(STIRAP) J4]. In this scheme, two atoms get combined to form a quasi-bound 
molecule in a higher electronic state via a FFR. The molecules formed get subsequently 
transferred to the ground state by two color laser pulses. In this situation, the total produc- 
tion rate is limited by the FFR conversion rate, the transfer efficiency of STIRAP being very 
high. Developing an understanding of the FFR phase of the ultracold molecule formation 
process thus remains important. 

The first reported experiment [s| swept the external magnetic field to move the unbound 
low temperature atoms through the FFR region. According to the experiments, the molec- 
ular conversion rate saturates at a very slow magnetic sweep, the so-called adiabatic re- 
gion. The adiabatic magnetic sweep is used for the ground state molecular production 
experiments ^] in order to maximize the number of molecules at a given temperature. This 
paper focuses on molecular production at the adiabatic region. The Landau-Zener(LZ) 
model has been used for analyzing the experimental results 5|-l8| . Detailed discussions of the 
LZ model and of its extensions are given in Ref.j^]. 

At any rate, Hodby et al.jf|, developed a semi-classical Monte Carlo simulation method, 
called the stochastic phase space sampling(SPSS) method to estimate the relation between 
phase-space density and molecular conversion rate by using a phase-space criterion depen- 
dent on a single parameter fitted to the experimental data at one temperature. No explicit 
reference to the LZ-type transition mechanism is made in this approach. Nevertheless the 
SPSS method has been applied to the 4O K- 4O K0, 85 Rb- 87 Rb 10], and 40 K- 87 Rb 



11 



12| cases, 



and the results show good agreement with the corresponding experimental data down to the 
lowest temperature realized in the lab. We note that recently an experiment on 40 K- 87 Rb 
molecular production has been carried out in a mixed Bose- Fermi system. The result of the 
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SPSS model, however, does not agree well with the experimental formation rate of molecules 
in the region of quantum degenerate temperatures)"" 



On the other hand, Williams et al. 



131 ] developed a theory of the FFR conversion rate 



based on a coupled Boltzmann-equation treatment that includes atom-molecule interactions. 
S. Watabe and T. NikunifjjJ extended the theory of Williams et al. and proposed the so- 
called (chemical) equilibrium theory (ChET for short in this paper). The ChET has been 
applied to noninteracting ideal atomic gases trapped by a harmonic potential, resulting in 
agreement with the 85 Rb- 87 Rb experiment down to the lowest temperature realized in the 



lab. And they commented in Ref. 



14|] on a possible difference between the SPSS model and 



the ChET, at temperatures below the BEC-critical temperature T c , a vital region that has 
not been experimentally examined in detail thus far. Let us note in passing that no explicit 
reference to the LZ-type transition mechanism is made in ChET either. 

As regards the SPSS model, no systematic assessment of its basic assumptions has been 
made. There has been no study of the meaning of the phase-space criterion adopted, nor of 
possible alternative criteria. Here we propose a few new phase-space criteria for comparison 
and analyze the temperature dependence of molecular conversion rate for each criterion. 
In doing so, we present an overview of the statistical evaluation of the ultracold molecular 
conversion rate through comparison between the SPSS and the ChET models. Furthermore, 
the molecular conversion rate below T c is estimated since its experimental and theoretical 
exploration is much desired. To understand the dependence on thermodynamic distributions 
of the atoms and on the phase-space criterion adopted, we choose the number-balanced K2 
system in a 3-dimensional harmonic oscillator as a prototype and analyze the features of 
SPSS in phase space. The results are compared with the ChET results. Through these 
comparisons, we clarify the distinct temperature dependence of the molecular conversion 
rate between the SPSS and ChET models qualitatively. Number-imbalanced systems are 
also considered for other effects such as the gravitational sag as well as the gap between 
two different values of T c . Incidentally, the ChET model with interacting atoms has been 
developed in Ref. 15J. This current paper, however, focuses on noninteracting systems only 
to understand the overall behavior. 

The paper is organized as follows. SecjTI] outlines the SPSS and ChET models and gives 
some additional background. Sec JIII Al discusses the results of the SPSS model for the 
prototypical K 2 case with four different phase-space criteria. Sec JIII Bl shows the results for 
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Rb2 and KRb cases, systems for which the temperature dependence of the conversion rate 
has been investigated experimentally. Sec JIII Cl compares the results of the SPSS model with 
the ChET. Sec JIVI concludes the paper. This paper uses T c = 0.94 fiwN 1 ^ 3 jks for the unit 
of temperature throughout, where h, uj, N, and ks correspond to the Planck constant, mean 
trapping frequency, number of bosonic atoms, and the Boltzmann constant respectively. 
And we define the "molecular conversion" rate as the number of formed molecules divided 
by the initial number of minority atoms, namely Xm = N molccu i c / N atomMnovity . 



II. THEORETICAL MODELS: SPSS AND CHET 



The theory of molecular formation in a cold mixture of two atomic species is based on 
the following effective Hamiltonian as a common starting point, namely 



H= / rfr[^S(r)^ a) i(r)+4(r)?{ M ^(r)] 

J 0=1,2 

+k J dr {<^r)^(r)^(r)+h.c.} 
+g A J dr ^(r)^(r)^ B (r)^(r) 

where i/)a and ipB are the atomic field operators for the two types of atoms, and 4>m is the 



molecular field operator |16|; 



"5? = -^ + 

and 

h 2 

H M = -TTf 1 + M r ) + e rcs 

2(m A + m B ) 

represent the single-particle Hamiltonian for atoms and molecules, respectively. Here the 
atomic masses and tub differ in general. U a {r) and Um{t) are the external potentials for 
the trap and gravity, respectively. The energy e res of the resonant molecule is tuned close to 
a collision threshold in experiments by varying an external magnetic field B. The constants 
k and qa pertain to the association/dissociation of an atomic pair into/from a molecule, 
and the atom-atom interaction, respectively. The LZ-type transition leads to k, but both 
SPSS and ChET bypass the use of the interaction term by employing physically-motivated 
assumptions. The atom-atom interaction term is dropped in what follows. In the mixture 
of two types of atomic gas, the interchannel transition for molecular formation is effective 
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when the magnetic field sweeps across the resonance energy. The experimentally observed 
and reported molecular conversion rate at a given temperature pertains to the saturated 
conversion rate that is measured at a reasonably slow magnetic sweep rate. Let us consider 
first the moment when the magnetic field is far away from the Fano-Feshbach resonance 
region so that the relevant classical Hamiltonian is that of an aggregate of independent 
atoms in the harmonic oscillator trapping potentials. The quantum distribution functions 
are given by 

fa = exp[(H a -fi a )/k B T]±l (1) 
where the Hamiltonian for a single atom of type a is 

H a = + l^y, (2) 

2m a 2 

where "+" corresponds to fermions, and "— " to thermal bosons. 

We extend the previously developed SPSS in order to treat the condensate bosons (BEC) 
using the truncated Wigner approximation to the phase-space distribution. The ChET on 
the other hand does not require a specific distribution function for BEC but rather, only the 
chemical potential. Let us employ the ground state wave function of the harmonic oscillator 
to represent a noninteracting BEC and apply the Wigner transformation using 

<j) (r) oc exp (-^^r 2 ) . (3) 



This yields 



2h 



P a (r,p) = ^gj (f+ f) ^ (r- f) e^' h dq 

K ex P ~— i~ r • ( 4 ) 



m a oj a h h 



This probability distribution does not 



depend explicitly on the number of atoms unlike the 



Thomas- Fermi(TF) wave function 10j. On the other hand, the number of BEC atoms in 



the harmonic oscillator trap is given analytically by 



rj~f \ 3 
Tr. 



N a , (5) 



where N®, and N a correspond respectively to the number of BEC atoms and the total 
number of atoms of type a at given temperature T. The probability of finding a condensate 
atom is then N®P a (r,p). 



Before going farther, note that the BEC transition temperature T c differs in general for 
the two types of atoms due to the mass difference. The TF (Thomas- Fermi) wave function 
for a repulsively interacting BEC tends to spread out more and get flatter near the center of 
the trap so that the distribution tends more toward smaller momentum values than in the 
noninteracting Eq. (j3J). To calculate the molecular conversion rate, Papp and Wiemanflo| 
represented the effect of BEC by assuming that if a pair of atoms find themselves within the 
TF radius, they then form a molecule regardless of their relative momentum. Thus Ref. [lo| 
does not generate a phase-space distribution. Instead, they use only the spatial extent of 
the TF distribution without setting any criterion on the momentum. 

A. SPSS: Stochastic Phase Space Sampling 



The SPSS is a method that samples candidate atoms for molecular formation by semi- 
classical Monte Carlo simulation. Ref. jfi] imposed a certain reasonable- looking phase-space 
constraint on the initial atom-pair distribution. The conversion rate is calculated under the 
following two assumptions. 

[1] Only the atomic pairs satisfying a given phase-space criterion can form a molecule. 

[2] Once a molecule is formed, it does not dissociate into atoms. Indeed, the lifetime of a 
molecule is known to be considerably longer than the time scale of an experiment 17 j. 




FIG. 1: (Color online) Conceptual diagram representing "Stochastic Phase Space Sampling 
(SPSS)". See steps [l]-[4] described in the text. 

A simulation is implemented in the following four steps, illustrated in Fig. [T] 

[1] Distribute atoms randomly in phase space in accordance with the thermodynamic 
equilibrium distribution at the initial temperature T, Eq.(TJ|. This temperature is 
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fixed, meaning that the process is regarded as isothermal. In the presence of a BEC, 
Eq.flS]), is used for the condensed atoms. 

[2] Search for atomic pairs that satisfy the phase-space criterion. Avoid double-counting 
by erasing the used atoms from the list of candidates. 

[3] Count the number of atomic pairs found in the random search. Regard it as the 
number of the formed molecules. 

[4] Return to Step [1], and repeat the steps until the statistical noise is reduced below 
the target noise level. 



Since the phase-space volume made available by a geometrical constraint plays a key role 
in SPSS, one of our goals is to assess the sensitivity of the method to the criterion adopted 
for pair formation. In what follows we define four independent phase-space criteria, the first 
one being the version initially employed by Ref-js] which yielded an extremely good fit to the 
experimental results for the Bose gas and for a two-component Fermi gas. As we will see in 
the discussion in Sec. 1111} all these separate criteria reproduce the temperature dependence 
of the experimental molecular conversion rate satisfactorily at thermal temperatures, but 
show some departures below the BEC critical temperature. 

Let us introduce the following four phase-space criteria; 
Criterion 1: 

|ArAp| < lq h, (6) 



Criterion 2: 



Criterion 3: 



Criterion 4: 



\Ar x Ap\ < jih, (7) 



(lAxA^IIAyA^IIA^I) 1 / 3 < lv h (8) 



\AxiApi\ < j s h (i = x,y,z) (9) 

where Ar corresponds to the spatial separation of an atom pair and Ap = m rm Av, where 
m rm equals the reduced mass, Av corresponds to the relative velocity, and j q , ji, j a , and 
7„ are the pairing parameters which are assumed to be independent of temperature and 
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density. Criterion 1 is introduced on the basis of the number of accessible quantum states 
in the relative coordinate. Criterion 2 is introduced under the premise that the low partial 
waves are important for molecular formation. The use of h instead of h for Criterion 2 is thus 
intentional. Criterion 3 is similar to criterion 1, but represents the full phase-space volume 
occupied by the pair. Criterion 4 is introduced as a counterpart to criterion 3, actually a 
most familiar form in statistical mechanics. In any case, there is no a priori numerical or 
theoretical support for any one of these criteria in the present context of molecular formation. 

As for the accuracy targeted in the present paper, the simulation program generates 
probability distributions, as given above by Eq.((T]), by the rejection sampling, and then 
searches for all pairs that meet the phase-space criterion. The period of the random number 
generator is about 2.3 x 10 18 for this study. This is sufficient for dealing with 10 6 atoms. 
The search ends when no remaining atomic pair satisfies the criterion. These searches are 
repeated a sufficient number of times, and the results are averaged. The error is of the order 
of the inverse of the square root of the number of atoms times the number of iterations. 
Throughout our calculations this error is reduced to within 1 %. We shall return to the 
discussion of these criteria in Sec. IIHI 



B. ChET: Equilibrium Theory 



The ChET was developed by E.Williams et al. 13J, and extended further by S. Watabe et 
al. [jjj]. It is based on the result of the coupled atom-molecule Boltzmann equation approach. 
In the ChET, one solves simultaneous equations consistently with equilibrium conditions to 
obtain the molecular conversion rate. The first equation is the conservation of the total 
number of atoms, 

N tot = N A + N B + 2N M (10) 
and the second equation states the constancy of the number of each atomic species, hence 

N b + Nm 

where Na, Nb, and Nm denotes the number of the majority atoms, of the minority atoms, 
and of the molecules, respectively. And the ratio a is defined by the initial ratio Nb/Na, 
where we consider Na = Nb , where Nm = thus 1 ^ a. There are two equilibrium 
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conditions, one representing chemical equilibrium, 

(J>A + fJ>B = fJ>M + 8 (12) 

where \i denotes the chemical potential of each component, and S denotes the detuning, 
namely the energy difference between dissociated atomic state and the molecular bound 
state. The other condition represents thermal equilibrium, T4 = T B = T M , where T denotes 
the temperature of each component. The population of each component is a function of ji 
and T. ChET uses two assumptions. One is that the magnetic sweep process is adiabatic, 
thus the total entropy of the system is conserved. And the other one is that molecular 
production halts at 5 =0, due to the conservation of momentum and energy. So the method 
first calculates the total entropy at 5 — > 00 at a given T and a, and then traces the adiabatic 
state until 5 =0. Thus the molecular conversion rate is given by equating the total entropy 
at 5 —> 00 and that at 5 =0, 

Sa(Too, ^oo,a) + Sb(Too, ^oo,b) = Sa(Tq, Ujq,a) + Sb(T , u,q,b) + Sm(Tq, /io,A + ^o,b), (13) 
where S denotes the entropy of each component. 



III. DISCUSSIONS 



Features of the molecular conversion rate are now explored for three statistically distinct 
types of atom-atom pair, namely Fermi-Fermi, Bose-Bose, and Bose-Fermi. To this end, 
we employ K 2 as a prototype for our analysis, considering 39 K for the boson, and 40 K for 
the fermion. The K2 system simplifies the situation on two accounts. Firstly, there is 
negligible gravitational sag thanks to the nearly equal masses of the isotopes. Secondly, we 



may presume that there are no more than two characteristic temperatures in the 



^2 system, 



namely the Fermi temperature Tp and the Bose-Einstein condensation point T c [18|]. It is also 



worthwhi 
40 K- 87 Rb 



ejpointing out that in a theoretical proposal as well as in a recent experiment with 



19j . the authors of the experimental study employed an optical dipole trap with 



the so-called "magic frequency" , which eliminates the gravitational sag. For this reason, the 
molecules 39 K- 40 K and 40 K- 87 Rb can be considered equivalent under an appropriate scaling 
of parameters [l9 ]. At any rate, the issues gravitational sag and the gap in T c are explored 
as supplementary items later in Subsection B. 
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Let us now compare the results of SPSS models and that of the ChET. To begin with, 
observe that the ChET is based on entropy conservation during an adiabatic magnetic sweep 
and on the condition that the statistical distribution of atomic pairs and that of molecules 
are always in equilibrium. The SPSS model presumes the conservation of the two-body 
local phase-space volume but no equilibrium. However, when atomic pairs are chosen the 
SPSS appears to be adiabatic since conservation of the phase-space volume amounts to 
adiabatic invariance and thus to conservation of entropy. But consistency with thermal 
equilibrium is not met. In the presence of an atomic BEC in the Bose-Bose and Bose- Fermi 
cases, ChET assumes the chemical potential of the bosonic atom component is 0, so that 
the conversion rate gets flat-lined. However, in the presence of two atomic BECs in the 
Bose-Bose case, ChET reveals certain temperature dependence of the conversion rate. For 
the sake of comparison, we determine the molecular conversion rate using the procedure 
developed in Ref. [b^] and review it in Sec. IIII CI 

A. Features with Prototypical K2 Molecules 

The following specific systems are considered. 

• Fermi-Fermi: 40 K(/ = §,m.y = — |)- 40 K(/ = = — ~) 

• Bose-Bose: 39 K- 39 K 

• Bose-Fermi: 39 K- 40 K 

Note that experimental data on K2 formation is known to us only for the Fermi-Fermi 



system 



6 . 



TABLE I: Experimental parameters from 



6] 



used in the present simulation. Nx is the number of 



K atoms, Wtrap the trap frequency, and 7's pertain to the criteria. 



N K ^trap 7Q 1L IV IS 

3x 10 4 2vrx(470, 470, 6.7) Hz 0.19 0.26 0.0085 0.050 



For the isotropic trap considered in this paper, we rescale the length and momentum 
according to the trap frequencies. The values of j q , 7/, j v , and 7 S are thus determined by 
the fitting to the experimentally available 40 K(/ = |,m/ = — |) 40 K(/ = |,m/ = — |) data, 
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and are used for the other Bose-Bose and Bose-Fermi cases as well. Table [J shows the values 
of various parameters employed in this section for K2. The accessible phase-space volume 
for the atomic pairs increases as 7 increases. The fitting parameters thus become smaller in 
the order of 7q > 75 > 7l/27t > jy. 

1. Temperature dependence of the conversion rate 

Stereotypical behavior of the conversion rate is summarily displayed as a function of tem- 
perature in Fig |2] for the four criteria altogether. The thermal limits look similar for the 




T/Tf T/Tc T/Tc 



FIG. 2: (Color online) SPSS molecular conversion rates for four criteria shown altogether. The 
red line corresponds to the result of Criterion 1, green to Criterion 2, blue to Criterion 3, and 
pink to Criterion 4. From left to right the panels correspond to 40 K- 40 K, 39 K- 39 K, and 40 K- 39 K, 
respectively. 

three types of systems. It is probable that this is due to the fact that the Hamiltonian and 
the geometrical constraints are all quadratic in r and p. And also the chemical potential 
plays little role in the thermal limits. Then the relevant quantities can all be expressed in 
terms of temperature-scaled distance and momentum, namely f = r/y/T and p = p/yT. 
The magnitude of the formation rate depends then on the value of 7/T which controls the 
phase-space volume. Specific temperature dependences become apparent when the chemical 
potential /x(T) plays a role, namely at T < T c for boson and at T <0.6Tp for fermion where 
the chemical potential changes sign. See Figj2j In the Bose-Bose case the conversion rate 
grows close to 100% as T — > 0, suggesting complete overlap of the two BECs in this limit 
which makes it easy to satisfy the imposed criterion. However, the slope of the molecular 
conversion rate is discontinuous at T c because of the BEC phase transition across this tem- 
perature. In the Fermi-Fermi case, the conversion rate rises monotonically as T lowers. The 
T=0 limit is finite, however. No discontinuity in the slope occurs at Tp since there is no 
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phase transition. (Note that the Fermi temperature Tp=1.93 T c is out of the range displayed 
in this figure.) The Bose-Fermi case shows a noticeable discontinuity in slope because of the 
BEC phase transition and of the subsequent diminishing overlap in phase space. 

2. Relevant single- and two-particle phase-space regions 

Let us begin this subsection with the following question: Should the four alternative 
criteria pick out similar regions of the phase space when jq, jl, jv, an d 75 are each made 
to fit to the observed molecular conversion rate? Let us look into the initial phase-space 
distribution of atoms that lead to molecular formation by comparing the results of the four 
criteria. 



Fermi-Fermi Bose-Bose Bose-Fermi 




Spatial Separation (10 5 a ) 



FIG. 3: (Color online) Master diagram for Q(Criterionl) showing phase-space distribution of 
molecule-forming pairs, p(Ar, Ap) Eq. Q3J Each consists of three panels corresponding to T > T c , 
T ~ T c , and T < T c . The triplet of distributions are arranged so that the left one is the Bose-Bose 
case, the middle one the Fermi-Fermi case, and the right one the Bose-Fermi case. The horizontal 
axis shows the relative spatial separation Ar, while the vertical axis shows the relative momentum 
m re d x Av. The brighter the color, the higher the density. j q =0.19 as in Table. HI 

To this end, we define the "two-particle" distribution function represented by the number 
of candidate pairs averaged over the implemented iterations, N iter , within a suitably cho- 
sen phase-space domain D(Ar, Ap) of a moderately small volume in the neighborhood of 
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Fermi-Fermi Bose-Bose Bose-Fermi 




Spatial Separation(10 5 a ) 



FIG. 4: (Color online) Master diagram for L(Criterion2) showing phase-space distribution of 
molecule-forming pairs, p(Ar, Ap) same as Fig. [3j 7 g =0.26 as in Table. HI 



Fermi-Fermi 



Bose-Bose 



Bose-Fermi 




u 1234012340123 

Spatial Separation (10 5 a n ) 



FIG. 5: (Color online) Master diagram for V(Criterion3) showing phase-space distribution of 
molecule-forming pairs, p(Ar, Ap) same as Fig. [3j 7 g =0.0085 as in Table. HI 



(Ar, Ap), 



p(Ar, Ap) = —— N[(A rj , A Pj ) e D(Ar, Ap)] 



(14) 



Figs. [3X 12 and [6] are our master diagrams showing p(Ar, Ap) comprehensively. Each 
panel consists of the three symmetry cases. Down each column, the temperature varies from 
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Fermi-Fermi Bose-Bose Bose-Fermi 




12340123401234 



Spatial Separation(10 5 a ) 

FIG. 6: (Color online) Master diagram for S(Criterion4) showing phase-space distribution of 
molecule-forming pairs, p(Ar, Ap) same as Fig. [3j 7 g =0.050 as in Table. HI 

T > T c to T < T c . And ao represents the Bohr radius. 

Next consider the physical implications of the way the candidate pairs are distributed. 
First, we consider how the candidate pair distributions differ, depending on the combination 
of atomic species. In the Fermi-Fermi case, the shape of the candidate pair distributions 
remains unchanged throughout these temperature regions. In the Bose-Bose case, we can 
see a change in the pair distribution below T c , namely a reduction in size toward the origin. 
This change comes from the condensation of the bose gases where their phase space volumes 
become very small. In the Bose-Fermi case, there are special features near zero temperature. 
And these features depend critically on the phase-space criterion. Above T c on the other 
hand, there are no marked differences in the candidate pair distributions for the three 
symmetry cases because they are all representable by the Maxwell-Boltzmann distribution. 
One major difference throughout the diagram is that Criteria 2,3, and 4 permit more widely 
spread distributions than does Criterion 1. Stated somewhat differently, the distributions 
simulated by Criteria 2, 3, and 4 are generally sparse with their maxima roughly in the 
regions cut away sharply by Criterion 1. The sharp cut-aways with respect to Criteria 2, 
3, and 4 do not introduce sharp edges in this Ar-Ap representation unlike Criterion 1. 
Incidentally, the converse is not necessarily true. Let us note that in the phase space of 
angular momentum and relative kinetic energy, both Criterion 1 and 2 show sharp edges. 
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The maximum value of the angular momentum corresponding to the edge is about 1.2H 
so that the s-wave contribution is dominant, but with nonnegligible p-wave contribution. 
In all the panels, the high concentration of most likely candidates appear along similar- 
looking hyperbolic-shaped arcs near the origin for the Bose-Bose and Fermi-Fermi cases if 
on different scales. 




Position R (10 5 a ) 

FIG. 7: (Color online) The upper two figures represent phase-space density distributions of indi- 
vidual components; (a) the fermion to the left and (b) the boson to the right at T=0.1T C . (c) 
shows the density of the molecule-forming pairs under Criterion 1. (d) shows the density of the 
fermion cut away by Criterion 1 when the boson is assumed to be concentrated entirely at the 
origin. These two resemble each other. 

On the other hand, the Bose-Fermi case exhibits a considerably different appearance. To 
aid in understanding, plotted in the two upper panels of Fig. [7] are the phase-space density 
distributions of the individual atoms in single-particle phase space. The one on the left 
represents the fermion, and the right one the boson which is condensed near the origin. As 
a result, the molecular formation is restricted to the narrow overlapping regions of these 
single-particle distributions. The lower two panels of the same figure mark the fermionic 
atoms that contribute to the molecular formation simulated by Criterion 1. The one on the 
left shows the distribution of the fermion atoms in the formed molecules which is extracted 
in the course of numerical simulations. The one on the right shows the same distribution 
but constructed by applying Criterion 1 to the single-particle distributions of the upper 
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two panels. We note that they look rather similar. The boson (BEC) distribution being 
concentrated near the origin, the criterion is naturally met either at short distances or at 
small momentum as in Fig. El (Bose- Fermi case at the lower temperature column.) 

There is another feature worth noting for SPSS. The likely distance between the atoms in 
the initial pair turns out to be much larger than the experimentally known final size of the 
molecule whose diameter is comparable to the scattering length at the final magnetic field. 
To demonstrate this point, Fig. [8] shows the reduced distribution obtained by summing over 
the momentum at T=0.1T C . The likely distance is considerably greater for Criteria 2, 3, and 
4 than for Criterion 1. This observation reflects certain traits of the SPSS approximation. 
The molecular formation takes place within the finite duration of the adiabatic sweep which 
ChET presumes is sufficiently long for thermalization. The SPSS prespecifies the phase- 
space volume that would be involved in the isothermal transport (which happens to be also 
"adiabatic" in SPSS because there is no change in the phase-space volume). The thermal 
distribution for each component of the final state can be evaluated at the temperature of the 
initial state, but then the system cannot be in chemical equilibrium. This is where the SPSS 
and ChET differ, yet the molecular conversion rate comes out surprisingly similar above 
T c . Let us note in passing that as the magnetic field is swept, the scattering length grows 
rapidly near the resonance. Across the unitarity limit, the scattering length grows without 
limit, causing even the widely separated pairs to interact. However, without dynamical 
calculations, it is not possible to assess the contribution of this resonant region. 



Fermi-Fermi Bose-Bose Bose-Fermi 




1 2 3 40 1 2 3 40 1 2 3 ' 

Spatial Separation(10 5 a ) Spatial Separation(10 5 a ) Spatial Separation(10 5 a ) 



FIG. 8: (Color online) Pair density distributions as functions of spatial separation Ar at T = 
0.1T C . From left to right, Fermi- Fermi, Bose-Bose, and Bose-Fermi. Each solid line pertains to the 
indicated Criterion. All the peaks representing the initial separation are significantly larger than 
the experimentally known final size of the molecule. 

Now consider very low temperatures below T c , as we explore the case T=0.1T C in what 



1(3 



follows. Figs. |3]|6] summarize the phase-space distributions under consideration. Fig. [8] plots 
pair density distributions as functions of the spatial separation p s (Ar) = J p(Ar, Ap)dAp. 
In the Fermi-Fermi case, the maximum appears near the origin, and then the number of 
candidate pairs decreases rapidly as the spatial and relative momentum separations increase 
away from the origin. This trend is independent of the applied criterion, but the height of the 
peak and the decline in distribution depend on the criterion. Because of the sparsely spread 
distribution as seen in Figs. HH6] away from the origin, the initial distance distributions 
of criteria 2, 3, and 4 are seen to decrease more gradually compared to criterion 1. In 
the Bose-Bose case, pair density distributions p s (Ar) show similar tendencies because two- 
particle distributions have the same trend, a consequence of the small size of the BEC for all 
the four criteria. In the Bose-Fermi case, a structure appears in addition to the peak near 
the origin. In Criterion 1, there are two peaks, one coming from atomic pairs with small 
spatial separation and high relative velocity, and the other one with small relative velocity 
and large spatial separation as noted earlier in relation to Fig. [71 In the other criteria 2, 
3, and 4, they also have double-peaked structure but the profile of these peaks strongly 
depends on the criterion. For example, the outer peak, somewhat broad and round, reaches 
the highest value for Criterion 3. This profile comes from the elliptic-shaped structure seen 
in the low right panel of Fig. |5j In any combinations of thermodynamic functions, all the 
peaks yield Ar much larger than the experimentally estimated final size of the molecule. 
SPSS merely enumerates the candidate pairs no subsequent dynamics is implemented. 

To summarize, we have explored the one- and two-particle distribution functions and 
have observed the following. 

• The two-particle distribution function behaves more or less the same under all the four 
criteria examined except for the Bose-Fermi combination at temperatures below T c . 

• For the Bose-Fermi system, this strong dependence on the imposed criterion is reflected 
in the molecular conversion rate. 

• Compared to the molecular size, the peak(s) in the two-particle distribution function 
corresponds to an atomic pair separation which exceeds the known size of the molecule. 
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B. Other Features 



Two aspects of the molecular conversion rate are now analyzed. One is the case where 
there exist two distinct BEC critical temperatures, and the other is the case of a differential 
gravitational sag that partially separates the two species. We begin with the former. As 
the temperature is lowered, the condensation of the component with higher T c introduces 
a decline in the molecular production since the other component still remains in a thermal 
distribution. In a real experiment, as was discussed in Ref. 10J with 85 Rb- 87 Rb, the second 
component fails to attain lower temperature due to technical difficulty and remains thermal. 
In such a case, the decline continues until T=0. See the left panel in Fig. |HJ The Fermi- 
Fermi case is not treated since no conspicuous features arise due to the absence of any phase 
transition across Tp. 




FIG. 9: (Color online) SPSS conversion rates for 85 Rb- 87 Rb using Criteria 1, 2, 3, and 4. The 
left panel assumes the same temperature for the two components. The right one simulates the 



situation where 85 Rb fails to get cooler at T=2.2T C 



To understand the implications of gravitational sag, it is instructive to see the phase-space 
density as a function of temperature for 40 K- 87 Rb. The BEC phase transition marks the point 
where the growing spatial gap between the two atomic species leads to almost nil conversion 
rate at lower temperatures. See Fig. [10] for the phase-space density representation of how 
the overlap evolves as T lowers. The gravitational sag can be made insignificant, however, 
by choosing an appropriate optical trap frequency referred to as the "magic" frequency in 
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FIG. 10: (Color online) The left panel shows the SPSS conversion rate for the K- 87 Rb case. The 
right column of panels shows the evolution of the phase-space density distributions for 40 K and 
87 Rb as functions of temperature. The upper one is for T > T c , and the lower one for T < T c . 
The rather large mass difference in 40 K- 87 Rb results in almost nil conversion rate toward T — > 
for Criterion 1. In this figure, the lines of the SPSS molecular conversion rate are scaled by 0.70 
to account for molecular losses that occur during rf association; here we followed the procedure 
described in Refs.flll. \v\. 

C. Comparison with ChET 



The ChET model is now applied to Fermi- Fermi, Bose-Bose, and Bose-Fermi cases with 
number ratio a set to either 1 or 2/15, and compared with the results of the SPSS model for 
criterion 1. The a = 1 case corresponds to the 
case corresponds to the Bose-Bose experiment 

100, 



ermi- Fermi experiment |6|] , and the a = 2/15 
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FIG. 11: (Color online) Temperature dependence of the molecular conversion rate for the Fermi- 
Fermi case. The red line corresponds to the result of ChET, and the green corresponds to the SPSS 
result, (a) the a=l result, (b) the a=2/15 result. In the figure for a=2/15, the characteristic 
temperature Tp corresponds to the Fermi temperature for the majority atom. 



The results for the Fermi-Fermi case are shown in Fig. [TTJ The ChET produces a con- 
version rate that increases with decreasing temperature. And at T — 0, it always reaches 
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100% and is independent of a. SPSS shows a similar trend, but indicates a saturation effect 
at near zero temperature as seen in Fig. [11] (a). Whether the saturation effect appears or 
not in the SPSS depends on a. 



S. Watabe et. al. 



20J also applied ChET with resonant interaction within the limita- 



tions of the mean field approximation. Their result indicates that the saturation effect is 
small. They concluded that the resonant interaction introduces a suppression of molecule 
conversion, which reduces the conversion rate somewhat from its maximal rate. 
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FIG. 12: (Color online) The same as Fig. 1111 except here for the Bose-Bose case. Upper panels cor- 
respond to S=0~ , lower to 5=0 + . For the o=2/15 case, T c corresponds to the critical temperature 
for the major bosonic component. In ChET, T c gets shifted and the plateau regions emerge. 

In the Bose-Bose case, both SPSS and ChET show similar behavior at thermal tem- 
peratures except for a subtle difference stemming from the relationship between the initial 
temperature T and the BEC transition temperature T c in the latter theory. See the encir- 
cled features in Fig. [T2l ChET assumes that the system follows the isentropic curve as the 
magnetic field is swept so that the final temperature differs from the initial one. This causes 
the difference in the number of atoms in the BEC as reflected in the shift of T c , in contrast 
to SPSS. The ratio r c =T c (t final )/T c (t initial ) can be simply evaluated near the two limiting 
cases. Setting T equal to the condensation temperature of the major atomic component at 
5 =0, and fis,o = Hm,o and applying an approximation to the chemical potential, we obtain 
from Eqfj3]at a ~1 
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These equations could allow a simple estimation of the shift in T c . Thus at a =1, r c =0.91 and 
at a =2/15, r c =0.99. Fig. [131 compares these approximate analytical results with numerical 
ones. This predicted region of a shift in T c , where the trend of the conversion rate suddenly 
changes, could provide an experimental test for the validity of ChET. 

1 

£0.9 
0.8 

0.2 0.4 0.6 0.8 1 

a 

FIG. 13: (Color online) Comparison between approximations and a numerical result for the ratio 
r c . At a ~1, both initial atomic clouds have a BEC component. However at a ~0, only the 
majority component has a BEC component. This causes differences in the ratio r c to appear as a 
function of a. 

Moreover, in ChET a plateau appears when the temperature lies between the two values 
of T c for a = 2/15 as displayed in Fig. [12j In ChET, the chemical potential of the majority 
atoms goes to 0, and thus the chemical potentials of the rest, i.e. that of the minority 
atoms and of molecules must become equal ( Eq. [T2l) . So in this situation, the conversion rate 
becomes independent of temperature while the rate attained depends on the trap frequency 
of the minority atoms and that of the molecules. Below the lower critical temperature, 
there are two possibilities; one is when the atomic BEC is more stable than the molecular 
BEC, and the other is the opposite situation. In the former case, the SPSS conversion rate 
decreases with decreasing temperature, and it goes to zero at T — > 0, and in the latter case 
it goes to 100% as T=0. And the results of ChET behave similarly. 

A marked difference between the two methods appears in the Bose-Fermi case as in 
Fig. [14j Consider a system with a =2/15 for which there are more bosons than fermions. 
The SPSS conversion rate begins to decrease at T c for both a = 1 and 2/15. However the 
rates of descent differ in two cases. For a = 1 the decline is steady, but for a = 2/15, the 
SPSS conversion rate is almost flat, decreasing very slowly with decreasing temperature; 
but it begins to drop suddenly at the temperature where the number of fermions exceeds 
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FIG. 14: (Color online) Temperature dependence of molecular conversion rate for the Bose- Fermi 
case, the same as in Fig. [TT1 In the ChET, shift of the T c and the plateau regions appear as in 



the number of thermal bosons. But in ChET the conversion rate gets flat-lined below T c . 
We emphasize that the temperature dependence of the molecular conversion rate comes out 
more or less the same for SPSS and ChET, but a noteworthy difference manifests itself in the 
Bose-Fermi case below T c , and it would be especially desirable to have more experimental 
tests available in this regime. 

IV. CONCLUSIONS 

This work has extended the SPSS analysis of the molecular conversion rate in a quadratic 
trap, including temperature ranges where no experimental explorations have been made to 
date. In particular, the temperature dependence of the conversion rate below T c has been 
considered at great length. This extension has introduced four distinctive geometrical con- 
straints for the pairing criterion in phase space that controls which atomic pairs contribute 
significantly to the molecular formation. Each constraint contains, apart from the geomet- 
rical information, a single parameter 7 which serves to fix the overall magnitude of the 
phase-space volume for molecular formation. Our study has examined the sensitivity of the 
molecular conversion rate to the geometry and to the magnitude of the parameter 7. A 
useful tool for analysis is the quantity we have denoted the "two-particle distribution func- 
tion", which is the probability distribution of the molecule-forming pairs in their relative 
coordinate |Af| = \?a — tb \ and relative momentum |Ap| defined similarly. Our investigation 
shows that once the parameter 7 is fitted to the experiment at T higher than T c , then the 
SPSS method reproduces the experimental results well in the known temperature ranges 
irrespective of the constraints (or criteria) used (Sec JIII Al . 



Fig.Hl 
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The benchmark system considered is K+K— > K2 in all the three symmetry combinations 
Fermi-Fermi, Bose-Bose, and Bose-Fermi, and the results of the calculations and analysis 
are the following. 

• The two-particle distribution function behaves more or less the same under all the four 
criteria examined except for the combination of Bose-Fermi at temperatures below T c . 

• For the Bose-Fermi system, this strong dependence on the imposed criterion manifests 
itself in the molecular conversion rate. 

• For all the three cases, the peak(s) in the two-particle distribution function corre- 
sponding to the separation of a pair of atoms exceeds the known size of the molecule. 

Thus, the Bose-Fermi case proved different because the single-particle distribution for the 
bosonic atoms and that for the fermionic atoms behave totally differently as they undergo 
condensation (Sec JIII A~2]) . Temperatures below T c remain a challenge for SPSS. 



The results of SPSS have been compared with those of ChET. This comparison exhibits 
a very similar temperature dependence of the conversion rate in Fermi- Fermi and Bose-Bose 
cases. One marginal difference is that in the Fermi-Fermi case the conversion rate of ChET 
goes to 100% at T — > 0, however in the SPSS model the saturation appears below 100%. 

In the Bose-Bose case, there are two limits as T — > 0, depending on the value of the BEC 
chemical potential, i.e. exothermic 5 < or endothermic 5 > 0. The molecular conversion 
rate tends either to 100% or to 0%. As a result, there is no substantial difference in the 
temperature dependence of the conversion rate between either theory. And in the Bose- 
Fermi case the conversion rate as a function of temperature has a different structure below 
T c . ChET yielded a flat-lined conversion rate whereas SPSS gives a slowly temperature- 
dependent conversion rate whose behavior is further dictated by the value of a. This paper 
does not consider effects of inelastic collision which may play an important role below T c . 
Further theoretical investigation is needed for deeper understanding of the dynamical behav- 
ior of the molecular conversion and an assessment between the theories and the experiments. 

This work has revealed what aspects of the atom distributions in phase space affect the 
conversion rate and how. The gained information may serve to give an insight when rigorous 
quantum transport theory is effected to visualize the phase-space evolution of the system 
during the magnetic sweep. After all, experimental effort for going below T c to examine the 
temperature dependence of the conversion rate appears seriously needed. 
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